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SUMMARY 


This  report  documents  two  ship  alrvake  turbulence  models  which  have  been 
developed  to  aid  in  the  analysis  and  simulation  of  the  shipboard  launch  and 
recovery  dynamics  of  conventional  and  VSTOL  aircraft.  The  two  computer  sub- 
routines represent  a large  deck  carrier  and  a FF  1052  frigate  air  wake.  The 
carrier  model  has  been  used  extensively  for  automatic  carrier  landing  simu- 
lations. The  frigate  alrwake  model  has  been  developed  more  recently  and  has 
not  been  widely  tested.  The  FF  1052  alrwake  model  has  been  Installed  as  part 
of  the  NAVAIRDEVCEN  VSTOL  Launch  and  Recovery  Program  described  in  reference 
(a).  The  model  has  been  analytically  validated  against  wind  tunnel  data,  but 
it  has  not  been  tested  in  a piloted  simulator. 

These  two  models  are  being  documented  in  the  hope  that  they  may  form  the 
basis  for  a standardized  method  of  simulating  turbulence  for  shipboard  launch 
and  recovery  analysis. 
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INTRODUCTION 

Two  turbulence  model  algorithms  have  been  developed  for  use  in  the  design 
and  analysis  of  conventional  and  vertical  take-off  aircraft.  These  are  an  attack 
carrier  model  (CVA)  and  a FF  1052  frigate  model.  The  CVA  model  has  been  used 
since  1971  in  the  simulation  of  manual  and  automatic  landings  aboard  large  deck 
carriers.  This  model  has  been  used  at  the  NASA  Ames  Research  Simulation  Facility 
in  piloted  simulations  of  the  F-14,  A-7,  and  EA-6B  aircraft.  References  (b)  and 
(c)  describe  the  data  and  equations  used  to  develop  the  model.  Free  air  turbu- 
lence, steady,  and  unsteady  burble  components  are  combined  to  produce  the  complete 
CVA  model.  To  make  the  algorithm  more  suitable  for  VSTOL  aircraft,  several  changes 
have  been  made  from  the  form  described  in  reference  (c).  These  are  as  follows: 

1.  The  ship  alrwake  velocity  components  are  programmed  to  exponentially 
decrease  in  magnitude  beyond  specified  bounds  of  altitude  and  lateral  position. 

2.  The  vertical  component  of  the  steady  air  wake  has  been  changed  to 
include  an  upwash  near  the  bow.  This  was  inferred  from  time  histories  of  AV8A 
launches  aboard  the  U.S.S.  Roosevelt.  (See  figure  1.) 

3.  An  option  to  provide  discrete  gust  sinusoids  in  place  of  the  ship  air 
wake  model  has  been  added.  These  inputs  follow  the  form  described  in  reference 

(f). 


The  model  is  limited  to  cases  where  the  wind  is  aligned  with  the  landing 
deck.  Wind  velocity  is  specified  by  setting  the  ship  speed  variable,  VS. 

The  second  model  has  been  developed  exclusively  for  VSTOL  aircraft  and 
is  intended  to  represent  wind  conditions  near  an  FF  1052  frigate. 

The  turbulence  algorithm  was  developed  from  data  derived  from  the  wind 
tunnel  test  of  a 1/50  scale  FF  1052  ship  model.  These  data  are  contained  in 
reference  (d).  The  structure  for  the  algorithm  was  based  on  work  performed  by 
the  Vought  Corporation  and  described  in  reference  (e) . However,  the  model 
described  herein  has  been  modified  from  the  original  format  as  a result  of 
additional  analysis  of  the  data  contained  in  reference  (d) . The  model  described 
herein  is  a somewhat  refined  version  of  the  Vought  model. 

This  model  was  programmed  on  the  NAVAIRDEVCEN  CDC  6600  computer.  The 
Fortran  listings  Included  in  this  report  represent  the  present  status  of  the 
two  models. 

The  subroutines  have  been  tested  in  batch  mode  runs  and  are  considered  to 
be  reasonably  accurate.  However,  the  author  takes  responsibility  for  any 
deficiencies  that  may  still  remain  in  the  code. 

The  FF  1052  air  wake  model  was  examined  by  computing  time  histories,  power 
spectra,  and  probability  distributions  and  then  comparing  these  with  the  wind 
tunnel  velocity  .time  histories. 

This  report  describes  sll  work  done  to  validate  and  adjust  the  FF1052  model 
to  agree  with  wind  tunnel  data.  Both  the  CVA  model  and  the  FF1052  model  are 
offered  for  inspection,  use,  and  criticism. 
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DESCRIPTION  OF  THE  CARRIER  AIR  WAKE  MODEL 

Appendix  A contains  the  Fortran  subroutine  listing  of  the  CVA  airwake  model. 
Appendices  B and  C provide  a definition  of  the  Fortran  variables  and  a flow  chart 
of  the  subroutine. 

Four  distinct  disturbances  are  combined  in  the  CVA  model  to  produce  a 
velocity  disturbance  in  each  of  the  three  aircraft  axis  directions.  These 
sources  Include  the  following: 

1.  Horizontal  and  vertical  velocity  deviations  which  result  from  vortex 
shedding  as  the  deck  pitches. 

2.  A steady  downwash  and  a horizontal  velocity  deficiency  which  vary  with 

distance  aft  of  the  ship. 

3.  Random  horizontal  and  vertical  airwake  velocity  components  which  vary 
in  RMS  amplitude  and  frequency  with  range  to  touchdown  and  wind  velocity. 

4.  Three-component  free  air  turbulence  which  is  independent  of  aircraft 
position.  The  model  is  representative  of  a large  deck  attack  carrier. 

A separate  discrete  gust  option  is  also  provided  which  generates  sinusoidal 
gusts  of  the  form  described  in  reference  (f).  This  option  may  be  selected  in 
place  of  the  airwake  model. 

In  addition  to  the  wind  model,  radar  position  error  is  included  which  is 
appropriate  for  simulation  of  the  SPN-42  automatic  landing  radar  system.  The 
simulated  position  error  increases  linearly  with  distance  from  the  ship. 

Each  random  component  of  the  model  is  generated  by  processing  the  output 
of  a pseudo-random  number  generator  vhiah  produces  a sequence  of  numbers  that 
approximate  a zero  mean,  unit  variance  Gaussian  probability  distribution.  The 
sequences  of  numbers  generated  are  uniquely  determined  by  the  integer  start  tip 
array  IRDST(I).  This  array  is,  in  turn,  determined  from  the  array  IRDSI 
(I,  Ionce)  depending  on  the  value  of  the  integer  flag  lonce.  lonce  may  take 
values  from  0 - 10.  This  technique  makes  it  possible  to  repeat  any  one  of 
10  turbulence  sequences  in  order  to  compare  different  aircraft  configurations. 

If  Ionce  **  0,  the  turbulence  generators  are  not  reset  after  each  run.  This 
produces  a long,  non-repeating  turbulence  sequence  to  gather  more  accurate 
statistics  on  a fixed  aircraft  configuration. 

The  calculation  sequence  of  the  subroutine  is  divided  into  trim  and  operate 
modes  depending  on  the  value  of  parameter  T Mode.  (I  Mode  = -1  in  trim,  and 
I Mode  ■ +1  in  operate.)  All  parameters  are  initialized  at  time  = 0 in  trim 
and  a new  set  of  turbulence  velocities  is  calculated  each  time  the  subroutine 
is  called  during  the  operate  mode. 

The  free  air  turbulence  component  of  the  model  is  taken  from  appendix  (b) 
of  reference  (h) . This  reference  specifies  power  spectral  gust  models  for  the 
three  velocity  components  as  follows: 
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* Aa 

Vo 


200/ vQ 3 


S 

Vo/lOO 


(ft/sec)2 
(rad/ sec) 


71.6/vo3 


(v^7Ioo  + x)(' 


-S 


Vo/100 


(ft/sec) 2 
(rad/sec) 


0 


The  corresponding  root-mean-square  gust  amplitudes  are: 
ug  * .76  M./sec  (2.5  ft/sec) 
wq  * .457  M./sec  (1.5  ft/sec) 
vq  * .728  M./sec  (2.39  ft/sec) 

The  frequency  parameters  of  the  free  air  turbulence  filters  are  determined  by 
the  trim  airspeed  vQ.  Since  the  algorithm  was  developed  for  conventional  flight 
simulation,  the  equations  do  not  permit  the  direct  substitution  of  zero  airspeed 
for  a hover  case.  For  low  trim  speeds,  the  parameter  is  arbitrarily  set  to  18.29 
M./sec  (60  ft/sec).  This  value  was  selected  as  a reasonable  estimate  of  the 
maximum  relative  wind  over  a moving  ship. 

The  turbulence  routine  options  are  determined  by  selecting  integer  flag 
parameters  as  either  0 or  1.  The  normal  values  of  these  flags  are  as  follows: 

Iwind  “ 1 

ILburb  **  1 

ILturb  ■ 1 

IBfrz  » 0 

Idisct  ■ 0 

Iubt  ■ 0 
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Ivbt 

Iwbt 


0 

0 


This  produces  a complete  ship  airwake  model  for  approach  simulation.  The 
complementary  values  produce  these  effects: 

Iwind  = 0 Return  all  zero  turbulence  values. 

ILburb  = 0 Select  alternate  downwash  function. 

ILturb  *=  0 Set  all  free  air  turbulence  components  to  zero. 

IBfrz  **  1 Fix  turbulence  range  variable  to  xd  = xfrz. 

Iubt  “ 1 Set  u component  of  airwake  burble  to  zero. 

Ivbt  = 1 Set  v component  of  airwake  burble  to  zero. 

Iwbt  = 1 Set  w component  of  airwake  burble  to  zero. 


Idisct  ■ 1 Select  discrete  gust  options  in  place  of  carrier  airwake  model. 


The  ship  dependent  portion  of  the  airwake  was  developed  primarily  from 
water  tunnel  tests  of  CVA  ship  models.  The  steady  components  consist  of  tabulated 
values  of  horizontal  and  vertical  velocity  functions  (Au/vs  and  Aw/vs)  versus 
range.  These  functions  are  shown  in  figure  2.  The  downwash  and  reduction  in 
relative  wind  is  determined  by  multiplying  the  tabulated  functions  by  ship 
speed  vs.  The  model  assumes  that  all  the  wind  is  generated  by  the  ship  and 
that  the  velocity  is  aligned  with  the  landing  deck. 


Periodic  wind  components  with  magnitudes  proportional  to  ship  pitch  are 
also  included.  These  sinusoidal  vertical  and  horizontal  components  propagate 
aft  from  the  ship  at  85  percent  of  the  ship  speed  and  decay  in  amplitude  with 

range. 

The  apparent  frequency  at  which  the  velocity  varies  is  determined  by  the 
aircraft  closing  speed  according  to  formula 


Magnitude  = cosine 


t + 


0.85 


♦ $ 


1.57 


where  o>p  = ship  frequency,  Av  = aircraft  closing  rate,  vs  3 ship  speed, 
xj  = distance  to  ship  eg,  and  if  = phase. 


The  magnitude  of  the  sinusoid  decreases  linearly  with  distance  aft  of  the 
ship.  The  burble  is  set  to  zero  outside  of  established  limits  fore  and  aft  of 
the  ship.  Additional  limits  on  magnitude  are  established  by  altitude  an  I lateral 
position  of  the  aircraft.  The  magnitude  is  assumed  to  decrease  exponentially  at 
distances  greater  than  100  feet  to  the  left  or  right  of  the  ship  and  at  altitudes 
greater  than  50  feet  above  the  deck.  Reference  (t)  indicates  that  the  burble  is 
largely  restricted  to  this  region. 
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Figure  2 - CVA  Ship  Burble  Steady  Wind  Ratios 
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The  random  number  sequences  used  by  the  turbulence  algorithm  are  generated 
by  calling  subroutine  Rand  once  for  each  element.  (See  reference  (a)  for  a list- 
ing of  this  routine.)  This  routine  generates  a random  sequence  by  taking  advantage 
of  the  roundoff  which  occurs  when  floating  point  numbers  are  represented  by  finite 
length  computer  words.  The  output  of  the  random  number  generators  are  passed 
through  a first-order  washout  filter  to  eliminate  low-frequency  components  and 
are  then  multiplied  by  a unit  sine  wave  which  shifts  the  signal  spectrum  to  the 
desired  frequency  range. 

The  random  components  of  the  burble  are  determined  by  passing  the  output  of 
the  random  number  generators  through  filters  having  the  following  form: 


ub2 (s) 
wb2(s) 


6(X)  • /2  • t(X)  • Input 
t (X)  • s + 1.0 

0.035  • VS  • /&. 66  • Input 
3.33  • s + 1.0 


The  variance  6(x)  and  the  time  constant  t(x)  of  the  u component  versus  range  are 
shown  on  figure  3. 

The  total  burble  is  computed  in  ship  axes  and  converted  first  to  inertial 
axes  and  then  to  aircraft  body  axes.  Free  air  turbulence  is  computed  directly 
in  body  axes  and  then  summed  with  the  burble  component  to  give  the  total  turbu- 
lence velocity  components. 

When  the  variable  Idisct  is  set  to  one,  discrete  sinusoidal  gusts  are 
calculated  in  inertial  axes.  Component  magnitudes  are  specified  by  the  variables 
Vmn,  Vme  and  Vmd,  and  the  frequencies  are  established  by  the  variables  Frth,  Frte, 
and  Frtd.  Samples  of  the  turbulence  components  are  shown  in  figures  4 and  5. 


DESCRIPTION  OF  THE  FF1052  TURBULENCE  MODEL 


A small  ship  airwake  model  has  been  constructed  using  tabulated  values  of 
mean  velocity  and  variance  in  velocity  as  described  in  reference  (d) . A computer 
listing  is  included  in  appendix  D,  and  a definition  of  program  variables  and  a 
subroutine  flow  chart  are  found  in  appendices  E and  F.  The  model  uses  linear 
superposition  to  combine  mean  velocities  with  random  velocities  to  yield  the 
total  turbulence  disturbance  at  the  aircraft  eg.  It  is  assumed  that  the  tur- 
bulence may  be  adequately  represented  as  having  a Gaussian  velocity  distribution. 

An  earlier  version  of  the  FF  1052  turbulence  model  described  in  reference 
(e)  used  second-order  filtering  of  random  inputs  and  cosine  extrapolation 
functions  to  represent  the  gust6.  However,  recent  analysis  indicates  that  some 
improvement  in  accuracy  and  a simultaneous  reduction  in  complexity  can  be 
achieved  by  using  first-order  random  number  filters  and  linear  and  exponential 
shaping  functions. 

Reference  (d)  summarizes  the  results  of  a wind  tunnel  test  performed  using 
a 1/50  scale  FF  1052  frigate  model.  The  turbulence  representation  was  developed 


Figure  3 - U Component  Burble  Time  Constant  and  Variance 
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Figure  4 - Sample  Burble  Components 
VO  - VS  - 18  M./Sec.  (60  Ft. /Sec.) 
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from  statistical  summaries  of  the  measured  turbulence  velocities  which  define 
the  mean  and  variance  of  the  three  velocity  components  recorded  by  an  array 
of  hot  wire  anenometer  probes  placed  downstream  of  the  ship  model.  The  mean 
and  variance  were  computed  from  a series  of  262  measurements  made  over  a 1.6- 
second  period  for  each  probe  location,  wind  speed,  and  wind  direction.  The 
turbulence  model  is  based  on  the  assumption  that  the  full-scale  turbulence  may 
be  accurately  represented  as  filtered  white  noise  with  the  proper  variance  and 
with  an  additive  mean  velocity  value  where  the  value  of  the  mean  and  the  varian 
vary  with  position  relative  to  the  ship  touchdown  bullseye. 

The  turbulence  model  assumes  a point  mass  aerodynamic  model  in  which  all 
aerodynamic  forces  and  moments  are  calculated  from  the  instantaneous  value  of 
airspeed,  angle  of  attack,  and  sideslip  at  the  aircraft  center  of  gravity. 

No  consideration  has  been  given  to  either  unsteady  aerodynamic  effects  or 
significant  variation  in  the  instantaneous  velocities  over  the  span  and  length 
of  the  aircraft. 

In  addition  to  the  basic  turbulence  model  structure,  provision  is  made 
for  extrapolating  the  model  beyond  those  regions  where  data  were  measured  in 
the  wind  tunnel.  Exponential  and  linear  extrapolation  functions  are  assumed 
for  those  regions  above,  behind,  and  to  either  side  of  the  region  where  data 
are  available.  The  air  is  assumed  free  of  turbulence  outside  of  the  extra- 
polation region.  Figure  6 illustrates  the  regions  where  the  turbulence  model 
applies  for  the  FF  1052. 

The  mean  and  random  turbulence  velocity  components  are  scaled  by  the  factor 
FI  which  is  formed  as  the  product  of  Fxl,  Fyl,  and  Fzl.  These  all  have  unit 
value  within  the  data  base  region  and  smoothly  approach  zero  outside  the  region. 
The  lateral  limits  of  the  air  wake  are  defined  by  Ymax  and  Ymin  which  locate  the 
bow  and  stern  of  the  ship.  The  test  data  base  extends  from  Y = -5.7  M.  (-18.75 
ft)  to  Y = 52  M.  (170  ft),  X = 0 M.  to  X - 76  M.  (2S0  ft),  and  Z = 4.4  M.  (14.58 
ft)  to  Z = 13.6  M.  (44.58  ft). 

All  distances  are  measured  relative  to  the  instantaneous  location  of  the 
touchdown  bullseye.  For  small  wind  over  deck  angles,  the  lateral  turbulence 
limits  are  defined  by  the  data  base  limits.  For  larger  heading  angles,  the 
tabulated  values  are  held  constant  and  extrapolated  to  Ymax  and  Ymin.  Beyond 
these  lateral  limits,  the  scale  factor  Fyl  follows  an  exponential  function 
decreasing  from  1.0  to  0.367  (e~l)  in  one  beamwidth.  Examination  of  the  tabulat 
data  shows  that  this  function  produces  a reasonable  match  in  the  region  to  the 
right  of  the  center  line  where  most  data  are  available.  Although  insufficient 
data  is  available  to  validate  the  assumption,  the  same  form  of  the  extrapolating 
function  ^e  Ho) /Beam!  used  to  extrapolate  vertically  beyond  the  upper  limi 
of  the  test  data.  This  function  is  consistent  with  the  shape  of  simple  boundary 
layer  flows.  In  the  streamwise  direction,  the  turbulence  limits  are  defined  by 
the  bow  for  zero  relative  heading  angles,  and  by  the  upstream  side  of  the  ship 
for  non-zero  relative  heading  angles.  The  extrapolation  limit  is  extended 
laterally  from  the  bow  and  stern.  All  turbulence  velocities  computed  for  points 
in  front  of  this  line  are  diminished  by  the  X extrapolation  function  Fxl  which 
decays  exponentially.  A linearly  scaled  extrapolation  function  is  used  for 
Fxl  aft  of  the  data  base  limit  of  76  M.  (250  feet). 
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The  amplitude  of  Fxl  decreases  linearly  from  a value  of  1.0  at  X = 76  M. 

(250  ft)  aft  to  0 at  X = 305  M.  (1000  ft)  aft  of  the  touchdown  point.  In  the 
regions  beyond  the  data  base,  the  nearest,  tabulated  value  is  used  with  the 
multiplicative  extrapolation  function  FI  applied  to  the  value  returned  from 
the  tables.  The  slope  of  Fxl  in  the  aft  region  was  suggested  from  plots  of 
variance  versus  downstream  range.  (See  figure  7.) 

The  random  components  of  the  turbulence  are  produced  by  filtering  the 
output  of  these  random  number  generators.  The  frequency  content  of  the  tur- 
bulence is  regulated  by  varying  the  interval  at  which  the  random  number  genera- 
tors are  called.  Increasing  the  calling  interval  decreases  the  frequency  of  the 
simulated  turbulence.  For  an  aircraft  simulation  interval  of  0.05  second,  a 
random  number  interval  of  0.3  to  0.35  second  was  found  to  produce  a reasonable 
turbulence  simulation.  The  random  number  sequence  is  linearly  interpolated  for 
time  values  between  succeeding  calls  to  produce  a smoothly  varying  number  sequence. 
This  sequence  provides  the  input  to  first-order  lags  which  effect  additional 
filtering  and  eliminate  high  frequency  variations. 

The  time  constants  of  the  filters  were  selected  by  examining  the  power 
spectra  of  samnle  test  data  points.  Significant  differences  in  the  frequencies 
of  the  three  velocity  components  were  uncovered. 

The  X component  was  found  to  have  a lower  effective  frequency  content  than 
the  Y and  Z components.  In  addition,  there  is  a definite  change  of  turbulence 
frequency  with  range.  It  was  determined  that  the  maximum  spectral  frequency 
decreased  with  range  downstream  from  the  landing  pad. 

Thus,  the  model  was  constructed  to  allow  the  three  filter  constants  to  vary 
linearly  from  X ■ 0 to  X = 58  M.  (190  ft).  This  appears  realistic  from  physical 
considerations.  As  the  turbulence  propagates  downstream,  rapid  velocity  fluctua- 
tions tend  to  be  diminished  by  shear  forces  whereas  more  slowly  varying  velocities 

tend  to  persist. 

The  random  velocity  components  are  combined  with  the  steady  mean  components 
to  yield  the  total  turbulence  velocities  in  wind  axes.  These  velocities  are 
then  transformed  into  inertial  axes  for  use  by  the  main  portion  of  the  aircraft 
simulation. 


ANALYSIS  OF  WIND  TUNNEL  DATA 

The  wind  tunnel  data  for  the  FF  105?  summarized  in  reference  (d)  were 
analyzed  to  determine  if  any  significant  trends  or  correlations  existed  among 
the  various  test  conditions.  An  attempt  was  made  to  reduce  the  mean  and 
variance  values  for  the  three  wind  speeds  tested  to  nondimensionalized  functions. 
To  achieve  this  end,  the  mean  and  variance  values  were  nondimensionalized  by  the 
freestream  velocity.  The  resulting  functions  were 

Vx  - V„,  Vy  Vz  6x  6y  and  _6z 

v ' v * V * v * v ’ v 

T OO  * oo  * oo  ¥ OO  v on  oo 

These  functions  were  plotted  versus  range  for  VM  •*  10,  18,  and  23  M./sec. 
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Figure  7 illustrates  these  functions.  It  can  be  seen  that  the  nondimens ionalized 
functions  for  the  three  velocities  agree  rather  well.  This  suggests  that  a gen- 
eralized model  could  be  used  which  would  be  valid  for  all  relative  wind  velocities. 
However,  the  present  model  is  restricted  to  the  three  freestream  velocities  used 
in  the  tunnel  tests,  and  a separate  set  of  nonlinear  tables  is  required  for  each 
of  these  three  velocities. 

Wind  tunnel  measurements  were  taken  for  model  roll  angles  of  -30,  0,  and 
+30  degrees.  These  data  were  investigated  to  determine  if  any  significant  trends 
existed  among  the  three  roll  angles.  Mean  and  variance  size  at  various  measure- 
ment stations  were  compared  for  the  three  roll  angles.  However,  no  definite 
functional  relationships  could  be  established.  Although  it  would  be  possible 
to  statically  interpolate  between  the  three  roll  angles  given  the  instantaneous 
ship  roll  angle,  there  is  no  evidence  to  suggest  that  this  model  would  be  more 
accurate  than  simply  using  the  values  for  zero  roll  angle.  Furthermore,  air 
operations  will  probably  not  be  attempted  with  roll  angles  near  ±30  degrees. 

Thus,  the  level  deck  wind  velocity  data  is  probably  the  most  applicable  of  the 
three  roll  angles. 

A dynamic  model  test  will  be  required  to  obtain  data  needed  to  develop  a 
turbulence  model  which  accounts  for  ship  motion.  It  is  likely  that  such  a model 
would  be  relatively  complex  in  order  to  account  for  the  periodic  variation  in  the 
turbulence  as  well  as  the  propagation  of  this  turbulence  downstream. 
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TABLE  I 

Selected  Wind  Tunnel  Test  Points  Used  in  Spectral  Analysis 
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TABLE  II 

Statistics  of  Frequency  Value  Corresponding  to  Maximum 
of  Spectra  Grouped  by  Velocity  and  Velocity  Component 


Freestream  Velocit 

y 

Free 

juency  - 

Rad/Sec . 

M./S. 

(ft/sec) 

Component 

Average 

Maximum 

Minimum 

RMS 

9.1 

30 

X 

60 

no 

30 

30 

9.1 

30 

Y 

76 

no 

10 

34 

9.1 

30 

Z 

90 

150 

30 

43 

16.46 

54 

X 

95 

190 

10 

52 

16.46 

54 

Y 

108 

260 

10 

99 

16.46 

54 

Z 

140 

290 

30 

90 

21.9 

72 

X 

90 

170 

30 

58 

21.9 

72 

Y 

130 

170 

70 

43 

21.9 

72 

Z 

290 

410 

150 

107 

TABLE  III 

1 


Statistics  of  Peak  Spectral  Velocity  Grouped  by  Velocity  Only 
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Because  the  number  of  samples  was  relatively  small,  the  power  spectra  were 
generally  very  erratic.  The  data  were  thus  smoothed  by  averaging  over  groups  of 
six  succeeding  frequency  points.  Figures  8 to  10  illustrate  samples  of  the 
smoothed  spectral  data.  The  sample  cases  shown  in  table  I were  examined  with 
regard  to  spectral  content. 

These  data  were  analyzed  to  determine  if  any  significant  trends  existed 
with  respect  to  the  spectral  distributions.  The  frequency  at  which  the  spectrum 
attained  its  maximum  value  was  computed  for  each  case  and  the  data  were  sorted 
according  to  freestream  velocity  and  velocity  component.  The  maximum,  minimum, 
average,  and  RMS  value  of  the  frequency  at  which  the  peak  of  the  spectrum  occurs 
is  tabulated  in  table  II. 

Grouping  the  data  by  freestream  velocity  only,  the  frequency  statistics 
can  be  summarized  as  shown  in  table  III. 

Figure  11  summarizes  the  frequency  trend  as  a function  of  test  condition 
velocity.  A least  squares  fit  of  the  three  gross  averaged  frequencies  yields 
the  relation 


max 

test 


3.837  + 2.2274  • V - rad/sec. 


Strouhal  scaling  based  on  model  length  is  used  to  obtain  the  equivalent  full- 
scale  frequency  as 


max 

full  scale 


0.07674  + 0.04548  • V - rad/sec. 


For  the  test  velocities,  this  becomes 


TABLE  IV 

Turbulence  Frequency  Versus  Airspeed 


Full  Scale 


1.440 

2.76 

3.53 


Ft . /Sec . 


30 

59 

76 


This  differs  from  the  second-order  filter  model  frequency  of  1.9  rad. /sec. 
at  V « 18  M./S.  (59  ft. /sec.)  as  described  in  reference  (d).  If  additional  wind 
tunnel  data  is  unavailable,  the  FF  1052  model  may  be  applied  to  larger  or  smaller 
ships  having  similar  geometry  or  the  same  ship  at  different  velocities  by  applying 
Strouhal  scaling,  where 


SN  - Strouhal  number 
f » Turbulence  frequency 
h **  Ship  beam 
V * Flow  velocity 
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Figure  8 - Measured  Gust  Spectra,  V = 9 M./Sec 
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Figure  xu  - Measured  Gust  Spectra,  V = 23  M./Sec. 
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The  turbulence  model  was  next  examined  to  see  how  the  properties  of  the 
simulated  turbulence  compared  with  the  measured  turbulence.  The  algorithm 
filter  parameters  were  varied  to  determine  if  an  improved  match  could  be  obtained 
between  the  simulated  and  test  data  turbulence.  The  basic  filter  equation  pro- 
posed by  reference  (d)  is 

6U.  y.  Z)  * / 4Cu)n3 

F(S)  - 

S2  + 2tiDnS  + wn2 


An  alternate  filter  equation  of  the  form 


F (S) 


6 • / 

S + a) 


was  examined  in  addition  to  the  original  formula.  Also,  a filter  of  the  form 


F (S) 


K(S  + a) 

S2  + 2(;tonS  + wn2 


was  investigated.  The  numerator  factor  as  well  as  the  frequency  and  damping  of 
the  denominator  were  varied.  The  baseline  values  suggested  in  reference  (d)  were 
u " 1.9  x V/59  and  C * 0.4. 


The  power  spectra  of  the  simulated  turbulence  were  computed  for  V * 23  M./S. 
Figure  12  illustrates  the  spectrum  for  the  X velocity  component.  It  is  apparent 
that  the  simulated  turbulence  spectrum  only  very  roughly  approximates  the  theoret- 
ical spectrum  which  corresponds  to  the  filter  parameters.  This  apparently  occurs 
due  to  the  limitations  of  the  digital  algorithm.  Further  comparison  of  figures 
10  and  12  show  that  the  measured  spectra  contains  more  power  at  high  frequencies 
than  does  the  simulated  spectrum.  This  agreed  with  the  previous  analysis 
indicating  that  the  peak  spectral  frequency  was  greater  than  the  natural  frequency 
of  the  algorithm  filter. 


The  algorithm  frequency  was  increased  in  an  attempt  to  match  the  peak 
frequency  derived  from  the  data.  Figure  13  shows  the  effect  of  increasing  the 
filter  frequency  to  / \ 

fa)Z  = 4.6  5.93  rad. /sec. 

The  spectrum  of  figure  13  has  a maximum  at  u 5 2.36  rad. /sec.  This  matches 
the  value  obtained  for  the  test  data  with  V = 23  M./S.  It  also  displays  an  in- 
creasing amplitude  versus  frequency  at  low  frequencies,  which  is  characteristic 
of  the  measured  data.  The  damping  parameter  of  the  filter  was  next  varied  to 
determine  its  effect  on  the  simulated  spectrum.  Figures  14  and  15  show  spectra 
with  <0  «*  5.9  and  £ * 0.1  and  0.7. 


Compared  to  the  case  with  5 ■ .4,  the  configuration  with  c c 0.7  had  its 
peak  amplitude  at  a lower  frequency  and  lower  amplitude  at  high  frequencies. 

The  c ■ 0.1  case  showed  a large  peak  amplitude  near  the  damped  natural  frequency 
and  a rapid  amplitude  roll  off  at  lower  and  higher  frequencies.  Thus,  neither 
configuration  appears  to  improve  the  correspondence  with  test  data  over  that 
achieved  with  the  original  value  of  C “ 0.4. 
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Figure  14  - Simulated  Turbulence,  0.7  Damping 
Z Velocity  Component 


Figure  15  - Simulated  Turbulence  Spectrum,  0.1  Damping, 
Z Velocity  Component 
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One  additional  filter  configuration  was  examined  in  which  a first-order 
numerator  term  was  introduced  with  a break  point  at  the  denominator  natural 
frequency.  The  resultant  filter  had  the  form 


F(S) 


K (S  + ui) 


S2  + 2 (0.4)  uS  + u2 


A representative  spectrum  for  this  configuration  is  shown  on  figure  16.  This 
configuration  shows  a more  pronounced  peak  and  a less  rapid  dropoff  at  high 
frequency  than  does  figure  12.  However,  it  does  not  appear  to  match  the  meas- 
ured data  as  well  as  the  configuration  shown  on  figure  13. 

Comparison  of  the  various  power  spectra  is  found  to  be  a less  than  conclusive 
method  of  establishing  the  proper  form  for  a turbulence  simulation.  The  measured 
data  was  limited  to  256  point  samples.  Because  of  the  extremely  random  nature 
of  the  turbulence,  the  short  data  sample  appears  inadequate  to  exactly  define 
the  power  spectra.  Even  the  artificially  generated  turbulence  is  found  to  have 
rather  irregular  spectra  although  the  equations  used  in  the  algorithm  should 
theoretically  produce  a smooth  spectra.  This  may  be  caused  by  either  the  finite 
samples  examined  or  by  the  approximations  of  the  digital  algorithm. 

Because  of  difficulty  in  defining  an  exact  mathematical  relation  for  the 
spectra  of  either  the  measured  turbulence  or  the  simulated  turbulence,  the 
modeling  effort  was  limited  to  matching  the  mean  value,  the  variance,  the 
frequency  at  which  the  power  spectra  attained  its  maximum  value,  and  the  overall 
slope  of  the  spectrum  above  and  below  the  peak  frequency.  Only  the  mean  and 
variance  can  be  matched  in  an  exact  manner  by  adjusting  the  mean  functions  and 
the  random  filter  gains.  The  remaining  characteristics  must  be  established  on 
a trial  basis. 


All  of  the  second-order  simulation  models  showed  a more  pronounced  peak 
and  a more  rapid  high  frequency  dropoff  than  did  the  test  data.  As  a result, 
it  was  concluded  that  a first-order  model  might  be  more  appropriate.  The 
principal  disadvantage  of  the  first-order  filter  is  that  the  output  velocities 
tend  to  be  more  jerky  and  less  smooth  than  those  produced  with  a second-order 
filter.  The  problem  was  recognized  as  being  caused  by  the  way  in  which  the 
random  number  sequences  were  generated.  When  a random  number  sequence  is 
generated  at  each  aircraft  integration  interval  (typically  0.05  second),  the 
bandwidth  of  the  random  noise  sequence  is  high.  The  magnitude  of  a purely 
random  sequence  is  essentially  uncorrelated  from  one  point  to  the  next.  In 
fact,  turbulent  velocity  variations  occur  over  finite  times,  and  real  measured 
data  shows  correlation  over  a very  short  time  duration. 

The  apparent  solution  to  this  problem  was  to  generate  the  random  number 
sequence  less  frequently  than  the  aircraft  equations  of  motion  were  integrated. 
By  interpolating  linearly  between  succeeding  random  numbers,  a smoothed  input 
was  generated  which  had  limited  bandwidth.  By  specifying  the  ratio  of  aircraft 
equation  update  interval  to  random  number  update  interval,  the  power  spectral 
content  of  the  simulated  turbulence  could  be  more  readily  adjusted  to  match  the 
test  data. 
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The  smoothed  random  number  sequence  can  be  further  processed  to  yield  a 
reasonable  simulation  of  the  measured  turbulence  data  by  passing  the  number 
sequence  through  a properly  designed  first-order  filter. 


Because  of  the  erratic  nature  of  the  measured  data,  it  was  decided  to 
examine  groups  of  points  at  various  downstream  ranges.  Individual  power  spectra 
were  normalized  by  dividing  by  their  peak  amplitudes.  Then,  these  spectra  were 
averaged  on  a component  by  component  basis.  The  data  selected  included  all 
cases  for  which  4>  = 0,  \p  = 50  deg,  V = 23  M./Sec.  The  same  process  was  performed 
for  all  data  for  which  <p  - 0,  ip  = 50  deg.,  and  V = 9 M./Sec. 


All  points  at  a given  downstream  range  were  lumped  together  regardless 
of  lateral  position  or  height.  Figures  17  to  20  show  the  effects  of  downstream 
station  and  freestream  velocity  on  the  character  of  the  normalized  spectra  for 
the  X,  Y,  and  2 velocity  components. 


All  four  figures  show  a more  pronounced  dropoff  of  spectral  amplitude  with 
frequency  for  the  X component  than  for  the  Y or  Z components.  Further,  figures 
17  and  18  show  an  overall  lower  frequency  break  point  than  do  figures  19  and  20. 
This  illustrates  the  fact  that  the  turbulence  frequency  decreases  as  the  distance 
downstream  from  the  ship  increases. 


A least  square  straight  line  curve  fit  was  first  applied  to  the  averaged 
test  spectra.  Then  the  break  frequency  of  the  first-order  filter  was  assigned 
a value  corresponding  to  the  point  at  which  the  straight  line  reached  a value 


of  0.5  times  its  value  at  zero  frequency, 
of 

i 1 


This  was  chosen  because  the  magnitude 


S + 

h) 


1 


reaches  0.5  at  S - u,  The  curve  fit  analysis  was  applied  to  the  cases  for 
V = 23  M./Sec.  (76  ft/Sec.). 


Based  on  the  assumption  that  Strouhal  scaling  could  be  applied,  the  filter 
frequency  was  made  proportional  to  freestream  velocity.  The  selected  frequency 
was  ratioed  to  the  value  obtained  for  the  V = 23  M./S  case  by  multiplying  by 
V/76.  The  selected  turbulence  algorithm  was  then  run  repeatedly  to  obtain  a 
good  estimate  of  the  spectral  properties  of  the  model.  These  spectra  were 
averaged  component  by  component  and  compared  with  the  test  spectra. 


The  filter  break  frequencies  selected  for  V «=  23  M./S.,  and  X * 0 were 
as  follows: 


TABLE  V 

Filter  Break  Frequencies  at  X = 0 
Component  Break  Frequency  (rad. /sec.) 


X 

Y 

Z 


8.17 

10.41 

10.0 


- 35  - 


: TdH-  rtfr'  _ 


A 


10  M./Scc 


Figure  17  - Normalized  Spectra  for  V * 10  M./Sec.,  X * 0 (Measured) 
X,  Y,  and  Z Components 


Figure  19  - Normalized  Spectra  for  V = 23  M./Sec.,  X = 0,  11  M.  (Measured), 
X,  Y,  and  Z Components 
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These  frequencies  are  varied  linearly  with  range  out  to  X «=  58  M.  (190  ft.), 
beyond  which  the  frequencies  are  held  constant  at  the  following  values: 


TABLE  VI 

Filter  Break  Frequencies  at  X c 53  M.  (190  ft.) 

Component  Break  Frequency  (rad. /sec.) 

X 1.58 

Y 3.0 

Z 2.88 

The  simulation  model  was  run  repeatedly  to  establish  a representative 
power  spectra.  These  averaged  spectra  were  then  plotted  against  the  averaged 
test  spectra.  Figure  21  shows  this  comparison  for  the  three  velocity  components 
with  X - 0 and  V " 23  M./Sec. 

Figure  21  shows  a reasonably  good  correspondence  between  the  scaled  test 
data  spectra  and  the  turbulence  simulation  spectra.  The  random  nature  of  the 
turbulence  makes  it  impossible  to  achieve  an  exact  spectral  match.  Therefore, 
the  comparison  between  the  test  data  and  the  simulation  data  is  not  totally  con- 
clusive. To  gain  additional  insight  into  the  nature  of  the  measured  turbulence, 
the  shape  of  the  velocity  probability  distribution  was  examined. 

Sample  test  data  were  selected  and  probability  distributions  were  computed 
for  groups  of  256  measured  sequential  velocity  values.  The  mean  was  subtracted 
from  all  values  to  produce  an  unbiased  sample.  Next,  the  maximum  and  minimum 
value  were  determined.  This  interval  was  divided  into  20  equal  divisions,  and 
the  256  values  were  sorted  to  determine  the  number  falling  into  each  interval. 
Finally,  these  individual  distributions  were  averaged  over  a number  of  test 
points  on  a component  by  component  basis. 

The  derivative  distribution  was  also  computed  by  determining  the  change  in 
velocity  from  one  time  interval  to  the  next,  and  by  sorting  these  values  according 

to  magnitude. 

A similar  analysis  was  performed  using  the  simulated  turbulence  algorithm  to 
generate  sequences  of  1024  points,  and  a comparison  was  made  between  the  simulated 
and  test  data  distributions. 

Figures  22  to  24  show  composite  velocity  distributions  for  the  three 
velocity  components  as  measured  in  the  wind  tunnel.  The  plotted  distributions 
appear  nearly  symmetric  and  approximately  Gaussian.  Only  the  X component  appears 
somewhat  non-Gausslan.  It  shows  a uniform  distribution  close  to  the  mean  with 
symmetric  falloff  to  either  side.  No  statistical  testing  has  been  performed  to 
test  the  hypothesis  that  the  distribution  is  Gaussian.  However,  it  appears  to 
be  a reasonable  assumption. 

Figures  25  to  27  show  the  equivalent  velocity  distributions  for  the  simulated 
turbulence.  These  plots  also  show  a good  approximation  to  a Gaussian  type  distri- 
bution. Thus,  it  may  be  concluded  that  the  simulated  and  test  data  velocity 
distribution  have  equivalent  forms. 


Figure  22  - Normalized  VZ  - Velocity  Distribution  - Composite  of  16  Test  Points 
V = 23  M./Sec.,  = 50°,  $ = 0,  X = 0,  27  M. 


d 


Figure  24  - Normalized  VX  - Velocity  Distribution  - Composite  Average  of  16  Test  Points 
VB  = 23  M./Sec.,  4 - 50°,  4 = 0,  X = 0,  27  M. 


500 


Figure  26  - Normalized  VY  - Distribution  - Composite  Average  - 4 x 1024  Points  - Simulated  Data 
V = 23  M./Sec.,  $ - 0,  = 50° 


Figure  27  - Normalized  VZ  - Distribution  - Composite  Average  -Ax  1024  Points  - Simulated  Data 
V - 23  M./Sec.,  <|>  « 0,  * « 50° 
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The  derivative  of  the  velocity  distributions  were  used  primarily  to  examine 
the  frequency  content  of  the  test  and  simulated  data.  The  variation  in  the 
velocity  and  the_variation  in  the  rate  of  change  of  velocity  were  computed. 

Then,  the  ratio  u = 6$/6V.  This  ratio  has  units  of  frequency  and  was  used  as 
a basis  for  matching  the  simulation  model  to  the  test  data. 

This  ratio  of  variations  was  computed  for  the  cases  used  to  calculate  the 
distributions  illustrated  in  figures  22  to  27. 


The  following  table  summarizes  the  comparison  between  test  and  simulated 
data  with  RNX  *»  RNY  = RNZ  = 5.0,  DTW  = .05  sec. 


TABLE  VII 

Test  Versus  Simulation  Equivalent  Frequency 


Component 


Equivalent  Frequency  - u 

Test  Data  Simulation  Model 


3.78  rad/sec. 
4.2  rad/sec. 
4.23  rad/sec. 


5.38  rad/sec. 
5.29  rad/sec. 
5.52  rad/sec. 


This  suggests  that  either  the  model  iteration  time  constant  or  the  filter 
time  constant  should  be  increased. 

A second  calculation  was  performed  to  compute  a mean  frequency  for  the 
test  and  simulated  data.  This  was  performed  using  the  spectral  data  in  the 
following  equation: 

N 

I S (io)  i * Au>  * ui 
. i = 1 


Average 


- rad/sec. 


E S(w)i  * Ait) 

i = 1 


This  is  the  eg  of  the  spectral  distribution. 


The  frequency  for  the  test  data  was  scaled  by  1/50  for  full-scale  comparison. 
This  averaged  frequency  parameter  showed  the  following  comparisons: 

TABLE  VIII 

Test  Versus  Simulation  Average  Frequency  (Averaged  Frequency  - rad. /sec.) 


Component 


Test  Data 


Simulated  Data 


4.7707 

4.4115 


4.7282 


- 48  - 
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The  fact  that  the  equivalent  frequency  of  all  three  components  of  the 
simulated  turbulence  is  nearly  identical  in  spite  of  differences  in  the  filter 
constants  suggests  t,hat  the  random  number  calling  interval  dominates  the  spectral 
distribution  of  the  simulated  turbulence.  Therefore,  the  program  scale  factors 
RNX,  RNY,  and  RNZ  were  adjusted  from  initial  values  of  5.0  to  RNX  ■=  7.12,  RNY  = 
6.3,  and  RNZ  ■ 6.52  based  on  the  ratios  of  the  equivalent  frequency  parameter. 

For  these  parameters,  the  equivalent  and  average  frequencies  compare  as  follows: 


TABLE  IX 

Test  Versus  Adjusted  Simulation  Frequency 


Test  Simulation 


Component 

0) 

uAv 

0) 

wAv 

VX 

3.78 

3.76 

4.51 

3.64 

VY 

4.2 

4.44 

4.63 

3.68 

VZ 

4.23 

4.54 

4.65 

3.79 

Given  the  accuracy  of  the  analysis  and  the  limited  sample  size  examined, 
this  appears  to  represent  a near  optimum  match.  Figure  28  illustrates  the 
resulting  spectra  for  a single  simulation  run. 


CORRELATION  ANALYSIS 


Certain  test  data  points  were  selected  for  correlation  analysis.  The  auto- 
correlation and  cross-correlation  functions  were  computed  to  determine  if  the 
gusts  were  purely  random  or  if  any  time  or  space  correlation  existed.  The 
equation  used  was 


As  expected,  the  autocorrelation  functions  were  found  to  have  unit  values 
for  zero  time  shift,  but  were  close  to  zero  for  all  non-zero  time  shift  values. 
This  means  that  no  significant  periodicity  exists  in  the  measured  data.  Most 
of  the  crosscorrelation  functions  were  found  to  be  essentially  zero  for  all 
time  shift  values.  The  condition  examined  was  X = 0,  Y = 0,  “ 50°,  $ e 0°, 

H * 4.2,  7.6  M.  There  was  no  space  correlation  between  velocities  at  H ■ 4.2  M. 
and  H « 7.6  M.  above  deck. 

The  only  significant  crosscorrelation  uncovered  was  between  the  X and  Y 
components  at  H = 4.2  M.  (14  ft.)  This  can  be  interpreted  as  a vorticity  in  the 
wind  velocity  about  the  Z axis  at  this  point.  This  corresponds  to  an  aircraft 


Figure  28  - Simulated  Spectra  With  Adjusted  Time  Constants 
X,  Y,  Z Components 
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TABLE  X 

Correlation  Summary  of  Test  Data  = 23  M./S.,  \p  = 50°,  <j>  = 0 

Point  1:  X = 0,  Y = 0,  Z = 4.2  M. 

Point  2:  X = 0,  Y = 0,  Z = 7 .6  M. 

AT  = 0.0061  second  model  scale 
= 0.31  second  full  scale 


Components 

AT-Shift  0 

1 

2 

Xl 

*1 

1.0 

0.288 

-0.057 

*1 

Yi 

1.0 

-0.08 

-0.1 

Zl 

Zl 

1.0 

0.064 

-0.014 

Yl 

-0.544 

-0.154 

0.094 

Yi 

Zl 

-0.05 

-0.05 

0.08 

*1 

Zl 

-0.1 

0.04 

0.018 

X2 

X2 

1.0 

0.24 

-0.01 

y2 

y2 

1.0 

-0.08 

-0.002 

z2 

z2 

0.998 

-0.058 

-0.206 

X2 

y2 

-0.303 

-0.129 

0.0409 

y2 

z2 

0.0974 

-0.081 

-0.095 

*1 

z2 

-0.1 

-0.001 

0.08 

*1 

X2 

0.237 

0.16 

-0.074 

Yi 

y2 

0.284 

0.0107 

-0.074 

Zl 

z2 

0.143 

0.0446 

0.0455 

Xl 

y2 

-0.15 

-0.076 

-0.076 

Yl 

z2 

0.075 

-0.174 

-0.0216 

*1 

z2 

-0.096 

0.121 

0.171 

51  - 


VPMPSMMMu'. 


1 





NADC-78182-60 


yaw  rate.  It  is  not  clear  whether  this  has  a significant  effect  on  aircraft 
motion  so  that  provision  for  including  this  effect  should  be  added  to  the  model. 
The  effect  does  not  appear  strong  enough  to  justify  additional  complexity  in 
the  model. 


SENSITIVITY  ANALYSIS  OF  THE  EFFECTS 

OF  THE  ASSUMED  FORM  OF  THE 

TURBULENCE  SPECTRUM 


Evaluation  of  the  test  data  and  the  various  model  algorithms  do  not  lead 
to  an  exact  and  definite  conclusion  regarding  the  values  of  the  model  parameters. 

A question  arises  as  to  whether  the  exact  turbulence  spectrum  must  be 
duplicated  by  the  model  in  order  to  yield  a useful  analytical  tool.  In  an 
attempt  to  answer  this  question,  the  closed-loop  pitch  response  to  horizontal 
gusts  was  computed  using  the  original  second-order  model  format.  RMS  pitch 
response  was  computed  for  a range  of  spectral  gust  frequency  model  parameters. 

A closed-loop  pitch  transfer  function  was  calculated  for  a representative  VSTOL 
(AV8A)  at  a 62  M./S.  (120  knot)  transition  condition.  The  pitch  loop  was  closed 
with  a simple  gain  pitch  SAS  and  pilot  model  incorporating  pitch  attitude  and 
pitch  rate  feedback  through  a delay. 


The  closed-loop  pitch  to  horizontal  gust  transfer  function  was  calculated 
and  the  RMS  gust  response  was  then  determined  using  the  equation: 


<5  e ^ 


i_  r i n uf 

!H  / D' 


|2  0 (jw) 

ugust 


where 


is  the  closed-loop  transfer  function 


and  ust  is  the  gust  power  spectrum. 


The  gust  power  spectrum  is  developed  by  squaring  the  magnitude  of  the  gust 
transfer  function  as  follows 


Gust  Velocity,  . _ 6 • 2 • /t;o)J 
Random  Input  ' S?  + 2coj  + o>2 


<J>  _(S) 

gust 


6^  A £(jj3 

(S2  + 2cwS  + u)2)  (S2  - 2cwS  + o)2) 


The  system  block  diagram  is  illustrated  by  figure  29. 


The  closed-loop  gust  transfer  function  then  becomes 

_6_,  . (3.28)  (S  + 8) (-O.OQ2832S2  + 0.0005604S  + 0.0003527) 

ugib'  “ (S5  + 7.046S1*  + 17.919983  + 27.827S2  + 16.692S  + 2.8202) 


1 
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NoS  = -8.057  S2  - 5.218  S - .823  Rad. /rad. 

D = S4  + 0.7382  S3  - 0.9896  S2  - 0.2859  S ♦ 0.006864 
Nug  = 3.28  X (-0.002832  S2  + 0.0005604  S ♦ 0.0003527) 

K0  = 0.42  rad/rad. 

Kq  = -0.49  rad/rad/sec. 

Figure  29  - Closed-Loop  Gust  Response  Transfer  Function 
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RMS  gust  responses  were  calculated  for  the  following  parameters 

a = 7.1  M./S.  (23.3  ft/sec.),  C = 0.4,  w = 1.4,  2,  3,  5,  and  7 rad/sec. 
ug 

The  results  follow: 


TABI.E  XI 

RMS  Pitch  Versus  Oust  Frequency 
06  - deg.  Gust  Frequency  - rad. /sec. 


3.16 

1.4 

2.95 

2.0 

2.43 

3.0 

1.76 

5.0 

1.44 

7.0 

For  the  data  case  analyzed  with  V = 23  M./S.,  the  original  model  uses 
a>x  = 1.9  x V/59  = 2.45  rad/sec.  Examination  of  the  data  suggests  using 
wx  *>  3 x V/18  = 3.86  rad/sec.  Interpolating  linearly  among  the  results  from 
the  sample  case,  we  get 


o8  - deg.  oi  - rad. /sec. 

2.72  2.45 

2.14  3.86 

This  corresponds  to  a 21-percent  difference  in  sigma  values  for  a 64-percent 
variation  in  frequency. 

Thus,  the  frequency  of  the  gust  model  does  significantly  affect  the  aircraft 
response  for  constant  RMS  gusts.  However,  a gust  model  which  predicts  gust 
response  to  within  10  percent  would  probably  be  a satisfactory  tool.  Thus,  the 
assumed  filter  algorithm  is  not  extremely  critical.  The  sensitivity  shown  by 
the  sample  indicate  that  it  would  be  sufficient  to  define  the  frequency  para- 
meter to  w’ithin  25  percent  of  its  true  value.  In  contrast,  the  aircraft 
response  varies  directly  with  RMS  gust  amplitude.  Thus,  it  is  most  important 
to  accurately  match  the  gust  RMS  amplitude. 


TIME  HISTORY  COMPARISON 


One  final  comparison  of  the  simulated  and  measured  turbulence  was  made  by 
plotting  the  velocity  time  histories  of  the  three  velocity  components.  The 
case  selected  was 

X = 0 V„  = 23  M./S. 

Y » 0 ^ - 50° 

Z - 4.2  M.  $ = 0 

- .4  - 


-y  * — 
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The  experimental  data  was  recorded  at  a time  interval  of  0.0061  second.  This 
value  was  scaled  up  by  a factor  of  50  to  yield  a plot  scale  interval  of  0.31 
second.  Simulated  data  points  were  calculated  at  a 0.05-second  time  interval. 
Because  of  the  measurement  interval,  the  bandwidth  of  the  experimental  data 
was  limited  to  10.13  rad. /sec.  Thus,  higher  frequency  oscillations  might  have 
been  present,  but  they  could  not  be  detected. 

The  simulated  and  test  velocity  sequences  were  compared  on  the  basis  of 
maximum,  minimum,  mean,  and  variance  values  as  well  as  a qualitative  comparison 
of  the  shapes  of  the  plotted  sequences.  Figures  30  to  32  illustrate  the  relation- 
ships for  the  three  velocity  components.  The  mean  and  variance  values  compare 
well  for  all  three  components.  The  maximum  and  minimum  values  seem  consistent 
in  all  cases  with  the  magnitudes  of  the  mean  and  variances. 

A qualitative  comparison  of  the  number  of  local  maxima  and  minima  in  a 
10-second  period  show  that  the  frequency  content  of  the  simulated  and  measured 
velocities  compare  satisfactorily. 


RESULTS  AND  CONCLUSIONS 

The  CVA  turbulence  model  contained  herein  remains  basically  unchanged  from 
the  version  described  in  reference  (c)  as  a result  of  analysis  performed  for  this 
study.  It  is  suitable  for  simulation  of  aft  approaches  and  forward  takeoffs  of 
conventional  and  VSTOL  aircraft. 

The  FF  1052  turbulence  model  is  designed  to  permit  VSTOL  takeoff  and  landings 
from  any  angle.  Wind  conditions  are  currently  limited  to  relative  wind  angles  of 
0,  30,  and  50  degrees  at  total  velocities  of  10,  18,  and  23  M./S.  (20,  35,  and  45 
knots) . 

After  considerable  iteration,  the  FF1052  turbulence  model  structure  was 
fixed  with  these  features: 

1.  Exponential  and  linear  extrapolation  of  tabulated  mean  and  variance 
velocity  values  beyond  the  region  of  measured  data. 

2.  First-order  exponential  filtering  of  random  number  inputs  using  inverse 
time  constants  that  decrease  linearly  with  range  downstream. 

3.  The  random  number  subroutine  used  to  generate  the  three  gust  velocities 
are  called  less  frequently  than  the  basic  simulation  loop  time,  and  the  resultant 
random  number  sequences  are  linearly  Interpolated. 

4.  The  random  number  update  interval  Is  selected  to  match  test  data  spectral 
distributions . 

5.  The  filter  time  constants  vary  inversely  with  freestream  steady  velocity. 

6.  The  frequency  content  of  the  three  velocity  components  display  different 
characteristics  with  the  horizontal  component  having  a lower  frequency  spectrum 
than  the  Y and  Z components. 


39.3  Ft/Sec 

23.3  Ft/Sec 


Component  Simulated  and  Measured  Turbulence 
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TIME  - SEC. 

Figure  31  - Y - Component  Simulated  and  Measured  Turbulence 
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7.  Although  some  cross-correlation  between  velocity  components  was  observed 
in  the  measured  data,  it  does  not  appear  significant  enough  to  include  in  the 
simulation  model. 

8.  For  45-knot  wind  and  the  aircraft  located  over  the  deck,  the  filter 
frequency  parameters  were  selected  with  the  following  values: 

Component  Frequency  - rad. /sec. 

X 8.17 

Y 10.41 

Z 10.00 

9.  These  frequencies  decrease  linearly  to  X = 58  M.  (190  ft.)  where  the 
filter  frequencies  become 

Component  Frequency  - rad. /sec. 

X 1.58 

Y 3.0 

Z 2.88 

10.  The  test  simulation  was  set  up  to  update  wind  velocity  intervals  at 
0.05  second  intervals.  The  random  number  sequences  were  updated  at  the  following 
multiple  interval  size: 


Component 

X 

Y 

Z 


Update  Interval  - Sec. 
0.05  x 7.12  = 0.356 
0.05  x 6.3  = 0.315 
0.05  x 6.52  - 0.326 


11.  Examination  of  the  experimental  data  shows  that  the  turbulence  veloci- 
ties reasonably  approximated  a Gaussian  distribution.  However,  the  distributions 
show  some  flattening  and  less  peaking  around  the  mean  than  do  ideal  theoretical 
distributions. 


12.  Analysis  indicates  that  it  is  more  important  to  exactly  match  the  mean 
and  variance  of  the  turbulence  velocities  than  to  exactly  match  the  spectral 
shape. 


RECOMMENDATIONS 

Because  of  differences  in  digital  computers  and  particularly  in  random 
number  generator  subroutines,  prospective  users  of  these  computer  subroutines 
should  follow  the  following  precautions: 


4 
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1.  Run  the  model  repeatedly  for  at  least  100  seconds  of  simulation  time 
at  a fixed  location  in  space. 

2.  Compare  the  mean  and  variance  of  the  simulated  data  with  that  of  the 
input  tables  corresponding  to  the  selected  point. 

3.  Compute  the  Fourier  transform  of  the  simulated  turbulence  and  compare 
these  with  the  spectral  plots  contained  in  this  report. 

4.  Adjust  the  filter  gains,  time  constants,  and  number  update  ratios 
RNX,  RNY,  RNZ  if  necessary. 

5.  The  FF  1052  model  should  be  tested  on  a moving  base  piloted  simulation 
to  get  pilot  confirmation  of  its  validity. 
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LIST  OF  SYMBOLS 

Attack  Aircraft  Carrier 
Ship  Heading  - Rad. 

Wind  Over  Deck  Relative  Angle  - Rad. 

Cross  Correlation  Function 
Laplace  Operator  - 1/Sec. 

Power  Spectral  Density 
Short  Take  Off 

Longitudinal  Wind  Velocity  - M./Sec. 

Longitudinal  Gust  Velocity  - M./Sec. 

Lateral  Wind  Velocity  - M./Sec. 

Lateral  Gust  Velocity  - M./Sec. 

Trim  Aircraft  Velocity  - M./Sec. 

East  Component  of  Natural  Wind  - M./Sec. 

North  Component  of  Natural  Wind  - M./Sec. 

Total  Wind  Over  Deck  - M./Sec. 

Ship  Speed  - M./Sec. 

X Component  of  Wind  Velocity  - M./Sec. 

Y Component  of  Wind  Velocity  - M./Sec. 

Z Component  of  Wind  Velocity  - M./Sec. 

Free  Stream  Wind  Velocity  - M./Sec. 

Lateral  Wind  Velocity  - M./Sec. 

Aircraft  Position  Aft  of  Ship  - M. 

Aircraft  Position  Left  or  Right  of  Ship  Centerline 
Aircraft  Altitude  Above  Touchdown  Point  - M. 
Aircraft  Angle  of  Attack  - Deg. 

Normalized  Horizontal  Velocity  Increment 

Normalized  Down  Wash  Velocity  Increment 

Aircraft  Pitch  Attitude  - Deg. 

Root  - Mean  Squared  Value 
Gust  Time  Constant  - Sec. 

Ship  Roll  Angle  - Deg. 


M 


r ^ 
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Horizontal  Gust  Power  Spectrum  - (f t /sec . )i.  2/ (rad/sec . ) 

Lateral  Gust  Power  Spectrum  - (f t/sec .) 2/ (rad/sec . ) 

Vertical  Cust  Power  Spectrum  - (f t/sec. )2/(rad/sec. ) 

Ship  Heading  Angle  - Deg. 


VO 


^o 
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APPENDIX  B 

Definition  of  Variables  Used  in  CVA  Program 

Common  Array 

Aircraft  Altitude  Above  Sea  Level  - Ft. 

Altitude  Boundary  on  Turbulence  - Ft. 

Common  Storage  Array 

Periodic  Amplitude  of  Sinusoidal  Airwake  Terms 
Time  Increment  - Sec. 

Nondimens ionalized  Time  Increment  of  u Burble  Component 

Nondimens ionalized  Time  Constant  of  W Burble  Component 

Time  Increment  - Second 

Extrapolation  Factor  of  Airwake  Velocity 

Extrapolation  Function  in  Y Direction 

Extrapolation  Function  in  Z Direction 

Shift  Frequency  - Cycle/Sec. 

Maximum  Gust  Frequency  - Rad. /Sec. 

Frequency  of  Vertical  Discrete  Gust  - Rad. /Sec. 
Frequency  of  East  Discrete  Gust  - Rad. /Sec. 

Frequency  of  North  Discrete  Gust  - Rad. /Sec. 

Common  Storage  Array 

Array  Containing  Turbulence  Gains 

Steady  Deck  Height  Above  Sea  Level  - Ft. 

Integer  Counter 
Integer  Common  Array 
Integer  Common  Array 
Constant  Range  Flag 
Discrete  Turbulence  Flag 
Downwash  Selection  Switch 
Free  Air  Turbulence  Switch 
Operating  Mode  Switch 

Turbulence  Initialization  Option  Switch 
Ship  Pitch  Flag 
Start  Up  Array 
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[ 


I 


4 


IRDST 

IuBT 

IVBT 

IWBT 

I wind 

IX 

OMEGP 

Phase 

PHIP 

PI 

PITCHM 

PITCIIMI 

RAV 

RD 

RDA 

RDX 

RDY 

RDYI 

RD3 

RD4 

RD5P 

RK10 

RN1 

SIN  WCT 

STHEIR 

SuBSW 

SuB2 

SwBSW 

S110 

S210 

TDTURB 

TIME 

TuB2 

Til  - T33 

u 


Start  Up  Array 

Switch  for  u Burble  Component 
Switch  for  V Burble  Component 
Switch  for  W Burble  Component 
Turbulence  Selection  Switch 
Random  Number  Startup  Parameter 
Ship  Pitch  Frequency  - Rad/Sec. 

Discrete  Turbulence  Phase  - Rad. 

Ship  Pitch  Phase  - Rad. 

3.14159 

Ship  Pitch  Magnitude  - Deg. 

Ship  Pitch  Magnitude  - Rad. 

Averaging  Parameter 

Filtered  Output  of  Random  Number  Generator 
Time  Constant  of  Random  Number  Filter  - 1/Sec. 

Raw  Output  of  Random  Number  Generator 

Intermediate  Random  Number  Value 

Intermediate  Random  Number  Value 

Intermediate  Turbulence  Variable 

Intermediate  Turbulence  Variable 

Intermediate  Turbulence  Variable 

Turbulence  Filter  Cain 

Radar  Noise  Gain 

Frequency  Shifting  Function 

Sine  of  Ship  Heading  Angle 

Steady  u Burble  Scale  Factor 

Random  u Burble  RMS  Magnitude 

Steady  w Burble  Scale  Factor 

Filter  Parameter 

Filter  Parameter 

Discrete  Turbulence  Switch  Time  - Sec. 

Simulation  Running  Time  - Sec. 

Random  u Burble  Time  Constant  - Sec. 

Euler  Angle  Conversion  Matrix  Relating  Aircraft  Body  Axes 
and  Inertial  Axes 

Output  of  Random  Number  Generator 
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| *• 

4 

i . 


uBSS 

UBSS1 

UBT 

uBl 

uBIA 

uBIK 

uBIS 

uBISS 

uBIST 

uB2 

uB2P 

uB2S 

uB2SS 

uB2ST 

uB2TA 

uG 

uGS 

uGSA 

uGSS 

uGST 

uTURB 

VBT 

VDBR 

VDTURB 

VEBR 

VEQIC 

VEQICP 

VETURB 

VG 

VGS 

VGSA 

VGSS 

VGST 

VG1 

VG2 

VG3 


Steady  u Component  of  Burble  - Ft. /Sec. 

Intermediate  Turbulence  Variable 
u-Body  Axis  Burble  Component  - Ft. /Sec. 

X-Ship  Axis  Periodic  Burble  - Ft. /Sec. 

Intermediate  Variable 
Intermediate  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 

Ship  X Axis  Random  Burble  Component  - Ft. /Sec. 

Intermediate  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 

Aircraft  X-Axis  Free  Air  Turbulence  Component  - Ft. /Sec. 

Statistical  Summary  Variable 

Statistical  Summary  Variable 

Statistical  Summary  Variable 

Statistical  Summary  Variable 

Total  Aircraft  X-Axis  Turbulence  - Ft. /Sec. 

Lateral  Burble  Velocity  - Ft. /Sec. 

Vertical  Burble  Component  - Ft. /Sec. 

Vertical  Turbulence  Velocity  - Ft. /Sec. 

East  Component  of  Burble  - Ft. /Sec. 

Aircraft  Trim  Velocity  - Knots 
Previous  Trim  Velocity  - Knots 
East  Component  of  Turbulence  - Ft. /Sec. 

Aircraft  Y-Axis  Component  of  Free  Air  Turbulence  - Ft. /Sec. 

Statistical  Summary  Variable 

Statistical  Summary  Variable 

Statistical  Summary  Variable 

Statistical  Summary  Variable 

Intermediate  Variable 

Intermediate  Variable 

Intermediate  Variable 
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VG4 

VMD 

VM£ 

VMN 

VNBR 

VNTURB 

VO 

VRW 

VS 

VSD 

VTURB 

VXBR 

VYBR 

WBSS 

WBSS1 

WBT 

WB1 

WB1A 

WB1K 

WB1S 

WB1SS 

WB1ST 

WB2 

WB2P 

WB2S 

WB2SS 

WB2ST 

WB2TA 

WC 

WG 

WGS 

WGSA 

WGSS 

WGST 

WTURB 

XCG 


Intermediate  Variable 

Vertical  Discrete  Gust  Magnitude  - Ft. /Sec. 

East  Discrete  Gust  Magnitude  - Ft. /Sec. 

North  Discrete  Oust  Magnitude  - Ft. /Sec. 

North  Burble  Component  - Ft. /Sec. 

North  Turbulence  Component  - Ft. /Sec. 

Trim  Velocity  - Ft. /Sec. 

Aircraft  Airspeed  - Ft. /Sec. 

Relative  Kind  Over  Deck  Speed  - Ft. /Sec. 

Burble  Propagation  Speed  - Ft. /Sec. 

Aircraft  Y Axis  Component  of  Total  Turbulence  - Ft. /Sec 
Ship  X Axis  Component  of  Burble  - Ft. /Sec. 

Ship  Y Axis  Component  of  Burble  - Ft. /Sec. 

Steady  Downwash  Velocity  - Ft. /Sec. 

Intermediate  Variable 

Aircraft  Z Axis  Component  of  Burble  - Ft. /Sec. 

Periodic  Component  of  Ship  Downwash  - Ft. /Sec. 
Statistical  Summary  Variable  - Ft. /Sec. 

Intermediate  Variable 

Statistical  Summary  Variable 

Statistical  Summary  Variable 

Statistical  Summary  Variable 

Vertical  Random  Burble  Component  - Ft. /Sec. 

Intermediate  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 
Frequency  Shift  Parameter  - Rad. /Sec. 

Aircraft  Z-Axis  Free  Air  Turbulence 
Statistical  Summary  Variable 
Statistical  Summary  Variable 
Statistical.  Summary  Variable 
Statistical  Summary  Variable 

Aircraft  Z-Axis  Turbulence  Velocity  - Ft. /Sec. 

Aircraft  X Position  Relative  to  Ship  eg  - Ft. 
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■t 


XD 

XFRZ 

X1S 

X2S 

YCG 

YR 

YRS 

YRSA 

YRSS 

YRST 

YR1 

YR2 

ZI1 

ZI111 

ZI112 

ZI2 

ZI3 

ZI4 

Z01 

ZOlll 

Z0112 

Z02 

ZR 

ZRS 

ZRSA 

ZRSS 

ZRST 

ZR1 

ZR2 


X-Position  Variable  Used  in  Table  Look  Up  - Ft. 
Fixed  Position  Used  in  Table  Look  Up  - Ft. 
Forward  Turbulence  Limit  - Ft. 

Forward  Turbulence  Limit  - Ft. 

Aircraft  Lateral  Position  Relative  to  Ship  - Ft. 
Lateral  Radar  Error  - Ft. 

Statistical  Summary  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 
Statistical  Summary  Variable 
Scale  Factor  on  Radar  Noise  Function 
Scale  Factor  on  Radar  Noise  Function 
Turbulence  Filter  Parameters 
Turbulence  Filter  Parameters 
Turbulence  Filter  Parameters 
Turbulence  Filter  Parameters 
Turbulence  Filter  Parameters 
Turbulence  Filter  Parameters 
Turbulence  Filter  Parameters 
Turbulence  Filter  Parameters 
Turbulence  Filter  Parameters 
Turbulence  Filter  Parameters 
Vertical  Radar  Noise  - Ft. 

Statistical  Summary  Parameter 
Statistical  Summary  Parameter 
Statistical  Summary  Parameter 
Statistical  Summary  Parameter 
Radar  Noise  Parameter 
Radar  Noise  Parameter 
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APPENDIX  C 

Flow  Chart  for  CVA  Model 


SINUSOIDAL  DISCRETE  GUST 
FRT  = AMAX1  (FRTN,  FRTE,  FRTD) 
PHASE  = FRT  * (TIME  - TDTURB) 

PHASE  > 6.28  Yes  


^ 1 

TIME  > TDTURB 


VNTURB*=J. 5 * VMN  * (1  - cos (FRTN  * (TIME  - TDTURB)) 
(37V-  RETURN 


NAUC-7M82-00 


COMPUTE  TURBULENCE  VELOCITIES 


COMPUTE  TOTAL  TURBULENCE  COMPONENTS 


UTURB  = UBT  + UG 

COMPUTE  AVERAGE  AND  RMS  VELOCITY 
PARAMETERS 
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APPENDIX  D 

Listing  of  FF  1052  Turbulence  Subroutine 
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APPENDIX  E 

Variable  Definitions  for  FF  1052  Program 


A 

ALT 

B 

Beam 

CPSIREL 

CPSISR 

CPSIW 

DT 

DTW 

DTX 

DTY 

DTZ 

EXPX 

EXPY 

EXPZ 

FOMEGX 

FOMEGY 

FX1 

FY1 

FZ1 

FI 

HTD 

IA 

IB 

IMODE 

IONCE 

IRDSTX 

IRDSTXI 

IRDSTY 

IRDSTYI 

IRDSTZ 


Common  Array 

Altitude  of  Aircraft  eg  Above  Sea  Level  - Ft. 

Common  Array 
Ship  Beam  - Ft. 

Cosine  of  Relative  Wind  Angle 
Cosine  of  Ship  Heading  Angle 
Cosine  of  Absolute  Wind  Direction 
Integration  Frame  Time  - Sec. 

Integration  Frame  Time  of  Turbulence  Subroutine  - Sec. 

Update  Time  Increment  of  X-Random  Number  Generator  - Sec. 
Update  Time  Increment  of  Y Component  Random  Number  Generator 
Update  Increment  of  Z Component  Random  Number  Generator 
X Component  Difference  Equation  Parameter 

Y Component  Difference  Equation  Parameter 
Z Component  Difference  Equation  Parameter 
X Frequency  Extrapolation  Function 

Y Frequency  Extrapolation  Function 
X Direction  Extrapolation  Function 

Y Direction  Extrapolation  Function 
Z Direction  Extrapolation  Function 
Total  Extrapolation  Function 

Height  of  Touchdown  Point  Above  Sea  Level  - Ft. 

Integer  Common  Array 
Integer  Common  Array 
Simulation  Mode  Switch 
Random  Number  Startup  Option  Switch 
Initial  Value  of  X - Random  Number  Sequence 
Current  Value  of  X - Random  Number  Sequence 
Initial  Value  of  Y - Random  Number  Sequence 
Current  Value  of  Y - Random  Number  Sequence 
Initial  Value  of  Z - Random  Number  Sequence 
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IRDSTZI 

Current  Value  of  Z - Random  Number  Sequence 

IWIND 

Turbulence  on/off  Flag 

IX 

Machine  Dependent  Random  Number  Generator  Parameter 

OMEGX 

X - Component  Filter  Constant 

OMEGY 

Y - Component  Filter  Constant 

OMEGZ 

Z - Component  Filter  Constant 

PI 

3.14159 

PSIREL 

Angle  of  Relative  Wind  With  Respect  to  Ship  - 

Rad. 

PS  IS 

Ship  Heading  - Deg. 

PSISR 

Ship  Heading  - Rad. 

PSIW 

Absolute  Wind  Direction  - Rad. 

RANDX 

Interpolated  Value  of  X - Component  Random  Number  Generator 

RANDXN 

Current  Output  of  X - Random  Number  Generator 

RANDXO 

Previous  Output  of  Y - Component  Random  Number 

Generator 

RANDY 

Interpolated  Value  of  Y Component  Random  Number 

Generator 

RANDYN 

Current  Output  of  Y .’landom  Number  Generator 

RANDY 0 

Previous  Output  of  Y Random  Number  Generator 

RANDZ 

Interpolated  Value  of  Z Component  Random  Number 

Generator 

RANDZN 

Current  Output  of  Z Random  Number  Generator 

RANDZO 

Previous  Output  of  Z Random  Number  Generator 

RLS 

Ship  Length  - Ft. 

RLTD 

Distance  From  Bow  to  Touchdown  Point  - Ft. 

RMODE 

Mode  Parameter 

RNX 

X Component  Time  Increment  Ratio 

RNY 

Y Component  Time  Increment  Ratio 

RNZ 

Z Component  Time  Increment  Ratio 

SIGVX 

X Velocity  Variance  - Ft. /Sec. 

SIGVY 

Y Velocity  Variance  - Ft. /Sec. 

SIGVZ 

Z Velocity  Variance  - Ft. /Sec. 

SIGX 

Tabulated  X Component  Variance  - Ft. /Sec. 

SIGY 

Tabulated  Y Component  Variance  - Ft. /Sec. 

SIGZ 

Tabulated  Z Component  Variance  - Ft. /Sec. 

SLOPEX 

X Extrapolation  Parameter 

SLOPE1 

X Boundary  Definition  Parameter 

SPSIREL 

Sine  of  Relative  Heading  Angle 

SPSISR 

Sine  of  Ship  Heading 
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SPSIW 

TCOUNT 

TCOUNTX 

TCOUNTY 

TCOUNTZ 

TFAC 

TIME 

VDTN 

VDTO 

VDTURB 

VDW 

VETN 

VETO 

VETURB 

VEW 

VEWBIC 

VNTN 

VNTO 

VNTURB 

VNW 

VNWBIC 

VS 

VW 

VW1 

VXBar 

VXT 

VXTP 

VXTURB 

VXWS 

VYBar 

VYT 

VYTP 

VYTURB 

VYWS 

VZBar 

VZT 

VZTP 


Sine  of  Wind  Heading 
Update  Time  Counter  - Sec. 

X Update  Time  - Sec. 

Y Update  Time  - Sec. 

Z Update  Time  - Sec. 

Turbulence  Update  Paramater 
Simulation  Time  - Sec. 

Current  Value  of  Vertical  Turbulence  Component  - Ft. /Sec. 
Previous  Value  of  Vertical  Turbulence  Component  - Ft. /Sec. 
Interpolated  Value  of  Vertical  Turbulence  - Ft. /Sec. 

Total  Vertical  Wind  Including  Free  Air  Wind  - Ft. /Sec. 

Current  Value  of  Computed  East  Turbulence  Component  - Ft. /Sec. 
Previous  Value  of  Vertical  Turbulence  Component  - Ft. /Sec. 
Interpolated  Value  of  East  Turbulence  - Ft. /Sec. 

Total  East  Wind  - Ft. /Sec. 

Initial  East  Wind  Component  - Ft. /Sec. 

Current  North  Turbulence  Update  - Ft. /Sec. 

Previous  North  Turbulence  Update  - Ft. /Sec. 

Interpolated  North  Turbulence  Component  - Ft. /Sec. 

North  Wind  Component  - Ft. /Sec. 

Initial  North  Wind  Component  - Ft. /Sec. 

Ship  Speed  - Ft. /Sec. 

Total  Steady  Wind  - Ft. /Sec. 

Velocity  Ratio 

Mean  X Turbulence  Velocity  - Ft. /Sec. 

Random  X Component  of  Turbulence  - Ft. /Sec. 

Previous  Value  of  VXT  - Ft. /Sec. 

Total  X Component  of  Turbulence  - Ft. /Sec. 

Tabulated  Mean  X Velocity  Component  - Ft. /Sec. 

Mean  Y Turbulence  Velocity  - Ft. /Sec. 

Random  Y Component  of  Turbulence  - Ft. /Sec. 

Previous  Value  of  VYT  - Ft. /Sec. 

Total  Y Component  of  Turbulence  - Ft. /Sec. 

Tabulated  Mean  Y Velocity  Component  - Ft. /Sec. 

Z Component  Mean  Turbulence  Velocity  - Ft. /Sec. 

Z Component  of  Random  Turbulence  - Ft. /Sec. 

Previous  Value  of  VZT  - Ft. /Sec. 
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VZTURB 
VZWS 
XE 
XMIN 
XMIN1 
XMIN  2 
XTaB 
XWIND 
YE 
YMAX 
YMAX1 
YMIN 
YMIN1 
YTaB 
YWIND 

ZTaB 


Total  Z Component  of  Turbulence  - Ft. /Sec. 

Tabulated  Mean  Z Turbulence  Velocity  - Ft. /Sec. 

North  - South  Distance  From  Aircraft  eg  to  Touchdown  Pt.  - Ft. 
Forward  Limit  on  Turbulence  Region  - Ft. 

Intermediate  Value  of  XMIN  - Ft. 

Intermediate  Value  of  XMIN  - Ft. 

X Table  Lookup  Parameter  - Ft. 

Aircraft  Position  Measured  Parallel  to  Relative  Wind  - Ft. 
Aircraft  East  - West  Position  Relative  to  Touchdown  Point  - Ft. 
Right  Limit  on  Turbulence  - Ft. 

Intermediate  Value  of  YMAX  - Ft. 

Left  Limit  on  Turbulence  - Ft. 

Intermediate  Value  of  YMIN  - Ft. 

Y Table  Lookup  Parameter 

Aircraft  Position  Relative  to  Touchdown  Point  Measured  Transverse 
to  Wind  Direction  - Ft. 

Z Table  Lookup  Parameter  - Ft. 
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1 


Ymin) 


FX1=  1 - (XTaB  - 250)  * Slope  X 


TIMl:  < TCOUNTX  ? 


RAND  X = RAND  X0+  (RAND  X N - RAND  X 0)  * 
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